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Abstract. The distribution of a time integral of geometric Brownian motion 
is not well understood. To price an Asian option and to obtain measures of 
its dependence on the parameters of time, strike price, and underlying mar- 
ket price, it is essential to have the distribution of time integral of geometric 
Brownian motion and it is also required to have a way to manipulate its dis- 
tribution. We present integral forms for key quantities in the price of Asian 
option and its derivatives (delta, gamma, theta, and vega). For example for 
any a > E [(At - o)+] = t - a + a 2 E [(a + At)- 1 exp(^^ - §)] , where 
A t = /* exp(B s - a/2) ds and M t = exp(B f - t/2). 
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1. Introduction 

The payoff of an Asian option depends on the (geometric or arithmetic) average 
of prices of a given risky asset over the pre-specified time interval. Under the Black- 
Scholes framework, one assumes that the price process {St, t > 0} of the risky asset 
follows 

dSt — fiStdt + uStdBt, So > 
where fi and a are given constants and {B t ,t > 0} is a standard one dimensional 
Brownian motion. In this setting, it is easy to understand the geometric average; 
if < <i < ti then in distribution 



y/S tl ■ S t2 = S exp 
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(aVt 2 + 3ti) • Af + (/i - +t 2 ) 



where Af is a standard normal random variable. On the other hand, the dis- 
tribution of an arithmetic average process is not well understood. A continuous 
version of the arithmetic average is a time integral of a price process. Using the 
Inverse Lapalce Transformation, Yor [16) proved many interesting identities re- 
lated to the distribution of geometric Brownian motion, which gives us deeper 
understanding of functions of geometric Brownian motion and useful information 
about their time integrals. More detailed research for the relation between the 
time integral and an Asian option was considered in [B]. Using the joint density of 
(J exp(B s )dW s , exp (.B t )) where B t ,W t are independent Brownian motions (given 
in [2]) the moment generating function of the time integral process was computed 
in [12] . The method of changing measures was considered to analyze the prop- 
erties of the time integral process. (See [7] [T3] [H] [IS]-) In [S] the very useful 
time reversing property is used to analyze the time integral process. Dufresne also 
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provided a certain form for the density function of the time integral of geometric 
Brownian motion. However the author pointed out the difficulties to use his for- 
mula in practice especially when the time integral is over the short time period due 
to the slow convergence rate. 

The payoff of an European style fixed strike Asian option is given by a function 
of the time integral of the price of the risky asset St 



where n is a fixed strike price and I is the pre-specified time interval with the length 
|/|. Under the risk neutral measure Q, we may set the price process St given by an 
SDE, dSt = crStdBt where a is a constant depending on the risky asset and Bt is a 
Q-brownian motion. Without loss of generality, we may also assume that the time 
interval I — [0, r] for some t > 0. It follows that the price of European style Asian 
option is given by 



i r 

- S cxp{er Bt ~ <r 2 t/2} dt - kY 
T Jo 



Since aB t = B a 2 t in law, we can rewrite the above quantity as follows 



(1.1) ^ rT E Q 



./ e B ^ 2 dt-^)+ 



Note that the initial price of the risky asset is not appeared in the time integral. To 
obtain the the price of Asian option and the derivatives with respect to the asset 
price we need to understand the following quantities 



(1.2) E[(A t -a)+], ^E[(A t -a)+], and ^E[(A t -a)+] 

where A t — J Q exp(B s ~s/2) ds and a > 0. In this paper we show that the quantities 
in (|1.2p can be expressed in terms of the expected values of functions of exponential 
Brownian motion. We believe that these expressions would provide the alternative 
approach to simulate the Asian option price and its greeks. The simulation results 
for Asian option price and its greeks were considered by several authors. (For 
example see [T], [3], g], [IT].) 

Other types of Asian options are also considered. When the strike price de- 
pends on the average price, it is called the floating-strike Asian option. (See for 
example [5] [5].) In [S], Henderson and Wojakowski show the very useful symme- 
tries between fixed-strike and floating strike Asian options. They showed that at 
the starting point of the averaging period there exists an equivalent relation be- 
tween the floating-strike Asian option and the fixed strike Asian option. However, 
once the averaging period has begun, the floating strike Asian option can not be 
re-expressed as a fixed strike option. 

In section 2 we discuss the relation between the time integral and the exponential 
Brownian motion. We summarize the result in [7] and present the key proposition. 
In section 3 we discuss the price of an Asian option and its derivatives (delta, 
gamma, theta, and vega). 
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2. The time integral of exponential Brownian motion. 



For /i G R>o let us denote Aj a time integral of an exponential Brownian 
motion with drift /i 

ri 



:= f exp[B a + (/i - l/2)s] ds. 
Jo 



When fi — we simply use A t without a superscript. For t > and y > we set 
M t := exp(B t - t/2) and 



^og(\- V -A t ). 
dt &y 2 ' 



y- 1 - 1/2 A t 

It is not hard to see that the process Rt satisfies the following SDE 

1 



(2.1) 



dR t = RtdB t - -Ri dt, R 



-y 



up to an explosion time t„ :— inf{t > : A t = 2/y}. By considering a stopping 
time t„ = inf{t : Rt < — n} and a Girsanov density process Rtl{t< Tn }i w e define a 
new measure Q under which 

r-l 

B t =B t 



y 



R s l{ s <r n } ds = B t - 21og(l - -A tATn ) 

is a standard Brownian motion. Let us define M = exp(B t —t/2) and A t = f Q M s ds. 
It follows that 

M t . v , 1 



(2.2) 
and 
(2.3) 



(1 



1 + |^tAT n 



1 - M« AT . 



RtAT n — ~ " 



M tAT „ 



y- 1 + iA tATre 

The following proposition is a simple modification of the proposition 3.3 in [7]. 

Proposition 2.1. If f(x, z) is a Borel measurable function and y > then 

2" 



E 



(2.4) 



/(M t ,A t );A t < 



/( 



y 



A t 



(l + fA t ) 2 ' 1 + 



] exp 



M 4 



2T 1 + \A t 



Proof. The proof of the proposition is essentially same as the proof of the proposi- 
tion 3.3 in [7], so we give a brief sketch. (For the detailed proof we refer [7].) For 
fixed n and y > it is not hard to see that 



SAt„ 



1 



*At„ 



R 2 S ds) 



exp(i? t + y)l{ Tn>t } = exp( / R s dB s 

Jo 

satisfies a Novikov condition. Thus we have 

E [/(M t) A t );r„ >t]=E Q [f{M u A t ) exp(-i? f - y)l TB>t ] . 

Since the event {t„ > t} is the same as {max s < t M S J (y _1 + 1/2 A s ) < n}, by letting 
ra -> oo the r.h.s. becomes Eq[f(M t , A t ) exp(-Rt - y)]. Use (|2^|) and we 
can rewrite the limit of the r.h.s. in terms of Mt,At and i?t. Since we are only 
interested in the quantity, we remove the tilde and get (|2.4ll . □ 
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Let us discuss certain of interesting choices of / in Proposition 12.11 Simple 
choices like a constant function or power functions allows us to have various relations 
between a simple martingale M t and its time integral. We consider the ones who 
directly related to the problem of pricing the European style Asian option and 
estimating the sensitivities. First we consider a constant function f(x,z) = 1. 
Using a constant function, it is easy to see that the l.h.s. of (|2.4| is a probability 
distribution of the time integral of a geometric brownian motion. By setting a = 2/y 
it follows that 

Lemma 2.2. For a positive a > 0, we have 
Pr[A t <a) = e - 2/a E 

For vel we have exp(t4/ 2 - v 2 t/2) = cxp(vB t - v 2 t/2) and thus 
(2.5) E\M% exp(j/t/2 - vH/2); A t < 2/y] = Pr[A { ^ ] < 2/y]. 

By setting f(x, z) = x" we get 

Lemma 2.3. For a real number v and a positive real a > 0, we have 



Pr[AY ] < a] 




Figure 2.1. The density function of A t . 

Let g t denote the probability density function for A t . It is known that the density 
function for At is continuous and positive. Figure 2.1 is the result of the simulation 
(10000 simulation for each a) using the identity given in the following Lemma [2^41 
The n-th moments of the time integral process At can be computed with simple 
computation. However it is known that A t has a heavy tail probability and thus 
knowing all integer moments does not give a probability density function. In Figure 
2.1 we can see the heavy tail probability. 

Lemma 2.4. For a real number a > 0, the probability density function g t for A t 
satisfies 

g t {a) = \pr[A t < a] ~ \pr[A^ ] < a}. 



exp 



M t 

a/2 + 1/2 A t 
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Proof. From Lemma 12.21 we have 



9t{a) = -rPr[A t < a] 
da 



\pr[A t < a] - \e- 2 / a E 



Mi 



■ exp 



(a/2 + 1/2 A t ) 2 \a/2 + l/2A t/ 

Since M t = exp(B t — t/2) we use M t as a Girsanov density function for the second 
term. Under the new measure B t — t/2 is a standard Brownian motion and thus 
the second term is equal to the following quantity: 



-e- 2 ' a E 
2 



1 



■ exp 



2exp(B t +t/2) 



a + A 



\a/2 + \/2A^Y 

Comparing with Lemma 12 . 31 we can see that the above quantity is the same as 
2/a 2 Pr[A i t 1) < a}. □ 

Furthermore we have 

Proposition 2.5. For any a > 0, b > and t > we have 

2M t \ 1 



Pr[M t <b,A t <a}= e" 2/a E 



exp 



a + A t 



; M t < 6(1 + -A t ) 



Proof. Consider f(x,z) = l{x < b} in equation (|2.4[) . By setting a = 2/y we get 
an integration form for the joint probability density function for M t and A t . □ 



The simulation result of the above proposition 12 . 51 for t — a is shown in the l.h.s. 
of the following figure and the joint distribution for t = 1 is shown in the r.h.s. 





Figure 2.2. Joint distribution of M t and A t . 

The first quantity of interests to have the derivatives of Asian options is the ex- 
pected value of the maximum function of the time integral subtracted by a constant. 
Using Proposition 12.11 we have the following result : 

Theorem 2.6. For any a > and v € R, we have 



E 



(A?) 



1 



l+e vt/2-„ 2 t/2 a 2„+2 E 



MX 



(a + A t ) 2v+1 



In particular if v = we have 
(2.7) E[(A t -a)+]=t- 



a 2 E 



(a + J 4 t )- 1 exp( 



2M t 

-i + At 



exp( 

_ 2, 
a' 



2M t 
a + A t 
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Proof. Let us set f(x, z) = x w (2/y-z). The l.h.s. of (J13J) becomes E[2/y-A t ; A t < 
2/y] = E[2/y - A t ] - E[(A t - 2/y)+]. Since 2/y - A t /(1 + y/2 A t ) = l/(2/y + A t ), 
the second result comes directly from (|2.4I) . Using the identity (|2.5p we can change 
the drift term of the exponential brownian motion to get the first result. □ 



Number of simulalio 



450 5Q0 



Figure 2.3. Expected value of max function : a = 0.4, £ = 0.5. 

The above figure is the simulation result of the second identity (|2.T[) in Theorem 
2.61 We set a = 0.4 and t = 0.5. The dotted line is obtained using the l.h.s. of 
the equation (12. 7p , that is we simulate At using Monte Carlo methods and take the 
maximum function. The solid line is from the r.h.s. of the equation (|2.7|) . Only 
with 500 simulation, we can observe that the convergence is faster for the solid line, 
which was expected because of the discarded simulations of the l.h.s. (when A t is 
smaller than a). 

Now let us consider the derivatives of E[(A t — a) + ] with respect to a. Since 



E[A t 



a; A t <a] = - / Pr[A t < u] du 
Jo 



we have that 



E[{A t -a)+} = E[A t -a}-E[A t -a;;A t <a] 

= t-a+ [ Pr[A t < u] du. 
Jo 



It follows that 

(2-8) ^E[(A t 
da 



= -1 + Pr[A t < a], and 



d 2 



a ) + ] = 9t(a). 
and Lemma [ 



The following theorem is the direct application of Lemma 

Theorem 2.7. The first and second derivatives of E[(A t — a) + ] with respect to a 
is given by following equations: 



-E[(A t o)+] = - 1 
da 

J^E\{A, - o) + ] . I 



-2/ajg 



exp 



2exp(B t - t/2) 



-2/a 



E 



exp 



a + A t 
2exp(B t - t/2) 



a + At 



AL 



M, 
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3. Sensitivities. 



Since At satisfies At+ S — At + MtA s for t, s > where A is an independent 
copy of A, pricing the fixed strike Asian option at time < t < the expiration date 
is essentially identical as the pricing the option at the beginning of the averaging 
period. 

One of important activities in Financial Market is managing the risk. One way to 
measure the risk in the option is estimating the ' Greek letters' such as delta, gamma, 
theta, etc. In this section we use identities obtained in Section 2 to get the ' Greek 
letters' of an European Style Asian option under the Black-Scholes framework; the 
price process of the risky asset follows under the risk-neutral measure 

dSt = aStdBt + rStdt, 

where Bt is a standard brownian motion under the risk-neutral measure, r is a 
constant interest rate and a is a constant volatility of given asset. The expression 
(jl.lj) is equal to the price of an European style Asian option with the expiration 
date r and the strike price k : 



Call — -^-e 

T<7 2 



r E 



e B t -t/2 dt _ O^KT^ 



-y 



(3.1) 

Theorem 3.1. The price of an Asian call at time t = is given by 

2M TCT 2 



Call = S e~ rT 



So 



S'2 
D 



-E 



(a + A Ta 2 ) cxp 



A T0 



whe 



So 



t is the expiration date, k is the strike price. 



Proof. The result comes directly from the combination of Theorem l2.6l and equation 
EH). □ 



The delta A is the rate of change of the price of the option with respect to the 
price of the underline asset. Thus we have 



(3.2) 



Call 



-E[(A a 2 T - a )+] 



So 



Sq Sq da 
Using Theorem 12.61 and Theorem 12.71 we have 

Theorem 3.2. The delta Aq of an Asian call option at time t = is given by 



A = e- rT + JL e -rr-2/a j a£ 



-E 



(a + A Ta 2 ) 1 exp 



2M TO 
i + A. 



exp 



2M T 



- A T(T 2 



So 



The gamma T is the rate of change of the delta with respect to the price of 
underlying asset. Thus the direct computation using equation (|3.2|) and the results 
in Theorem 12.61 Theorem 12.71 and Theorem 13.21 we have 



Theorem 3.3. The gamma Tq of an Asian call option at time t = is given by 

r o = " 7a ^ TT g Ta 2(- 



2 2 

a k r 



,a 2 KT, 



So 



where g t is a continuous probability density function of A t given in Lemma 
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Time lo expiralion 



Time to Expiralion 



Figure 3.1. Derivatives vs. the time to expiration: So = K = a = 1 

In Figure 3.1. we use Monte Carlo methods to get Greek letters A,r, 9 and v. 
We take the risk free rate r = and both the initial stock price and the strike price 
are equal to 1. We also set for simplicity £7 = 1. Notice that we plot greek letters 
vs. time to expiration. Since we computed greeks at the beginning of the averaging 
period, time to expiration is the same as the length of the averaging period. In the 
upper l.h.s. figure we can observe that as time to expiration increases, the value 
of Ao increases. The plot of Tq with respect to time to expiration is given in the 
upper r.h.s. figure. Since T is the first derivative of A with respect to the stock 
prices, it reflects our observation in the plot of Ao- The lower l.h.s. figure is the 
plot of @o and the lower r.h.s. figure is the plot of Vq. 

Since all parameters are constants, the call price satisfies the differential equation 



dCall 
dt 



■S 



dCall 1 2 B 2 Call 



It follows that 



6o + tSqAo + -o- 2 SoT = rCall 



where Oo is the rate of change of the price of the option at time t = with respect 
to time and Callo is the price of the call option at time 0. 



Theorem 3.4. The theta Go of an Asian call option at time t — is given by 
9 = rCallo - rS A - ^o- 2 S T 

where Callo, Aq and To are given in the previous theorems 3.1-3. 

Now let us discuss the vega v. The vega v is the rate of change of the price of 
the option with respect to the volatility of the underlying asset. The call price of 
European style Asian option is given in (|3.ip . Thus we have 
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dCall 



2 r> 11 i - 
■— Call H e 

cr cr 



■ CaZZ 



2 ^o . ,-- , 



cr a 

25Q _ rr 



< " E 



M a 2 T 



hi . a 2 KT 

— ■ A 2 > 



hi 



So 



M<j2 T - — ; A a 2 T < — — 

DO ^0 



It follows that 
dCall 



(3.3) 



dc 



2 2 9 
CWZ + 



1 -£ 



^cr2 r ; A a 2 T < 



S 



A a 2 r < 



So 



By setting v = 1 in equation (12. 5| in Lemma l2~2l we see that 



1 - E 
Therefore we have 



M a 2 T ; A a 2 T < 



So 



P 



A^ > 



a 2 KT 



Theorem 3.5. The vega vo of an asian call option at time t = is given by 

2k 



vo 



(T cr 



Ait> 



a 2 KT 
So 



-P 



A a 2 T > 



a 2 KT 



Last two quantities in the above theorem can be obtained using Lemma 12.31 
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